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Abstract. It is widely accepted that quasars and other active galactic nuclei (AGN) are powered 
by accretion of matter onto a central supermassive black hole. While numerical simulations have 
' demonstrated the importance of magnetic fields in generating the turbulence believed necessary 
for accretion, so far they have not produced the high mass accretion rates required to explain 
MH . the most powerful sources. We describe new global 3D simulations we are developing to assess 
' the importance of radiation and non-ideal MHD in generating magnetized outflows that can 
enhance the overall rates of angular momentum transport and mass accretion. 
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5/5 | 1. Introduction 

i Accretion onto a black hole is widely believed to power the central engine of high- 

energy astrophysical objects such as active galactic nuclei (AGN) and some X-ray bina- 
ries. Since the standard theory of astrophysical disk accretion was formulated over 30 
years ago (Shakura fc Sunyaev 1973 INovikov fc Thorne 1973|) . arguably the most impor- 



I ' tant advances in understanding these systems have come from computational modelling. 
. Numerical simulations demonstrate unequivocally that the magnetorotational instability 
(MRI) can produce magnetohydrodynamic (MHD) turbulence and the outward transport 
^vq , of angular momentum required for accretion to proceed (see IBalbus 20031 for a review) . 

Notwithstanding the important advances made using MHD simulations over the last 
ON . two decades, they have so far been unable to resolve two key outstanding issues: 1. How 
are the high mass accretion rates inferred in the most powerful sources achieved? and 2. 
, How are the outflows and jets observed across the mass spectrum of accreting sources 
produced? Our new work will test the hypothesis that these issues may be connected and 
mutually resolved by a generalized model for MHD disk accretion. 



2. MHD Disk Accretion Theory 

The analytical foundation for our new work is based on the generalized accretion disk 



model of Kuncic & Bicknell (2004) which is schematically illustrated in Fig. Q] The model 
includes contributions from both radial and vertical transport of angular momentum to 
the overall mass accretion rate: angular momentum is transported radially outwards by 
internal MHD stresses and vertically outwards by both mass outflows in a wind and MHD 
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Figure 1. A schematic illustration of the inner regions of an MHD accretion flow around a 
black hole. /3 is the usual plasma beta (the ratio of gas to magnetic pressures). The wind can 
drive mass loss from the disk, while the jet may be dominated by Poynting flux. 

stresses acting over the disk surface. This vertical transport of energy and momentum 
is not modelled in standard accretion disk theory. Our new simulations are designed to 
compare its effects with the radial turbulent transport of angular momentum and to test 



the prediction of a modified emission spectrum (see e.g. Kuncic & Bicknell 2007a,b) 



3. New MHD Accretion Simulations 

To implement our model, we are using and extending FLASbfj] ( |Fryxell, et al. 2000[ ), 
the public MHD code developed at the University of Chicago. FLASH implements adap- 
tive mesh refinement and solves the equations of MHD in conservative form thus explicitly 
conserving energy. 

To date, most simulations of the MRI in accretion disks have been conducted in a 
shearing box approximation but this approach introduced complications and pitfalls (see 
e.g. |Regev fc Umurhan 2008| . We will instead conduct global 3D simulations. 

By explicitly modelling a finite resistivity, we will explore its effect on the evolution 
of the magnetic field topology, particularly the self-consistent emergence of a significant 
z-component necessary to produce outflows and high mass accretion rates. 

We will include radiation in our simulations to directly compare our results to the ob- 
servational data. This will allow us to test whether the blackbody emission from the disk 
is modified by outflows as predicted. Radiative transfer is also be required to transport 
the internal energy dissipated in the disk plasma at the end of a turbulent cascade, i.e. to 
cool the disk. In addition, it may be that regions of the disk where radiation dominates 
may be thermally unstable ( |Shakura fc Sunyaev 1976] ), thus affecting the dynamics. 
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f FLASH is freely available at http://flash.uchicago.edu. 



